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ABSTRACT 



We revisit the possibility that we inhabit a compact multi-connected flat, or 
nearly- flat. Universe. Analysis of CO BE data has shown that, for such a case, the 
size of the fundamental domain must be a substantial fraction of the horizon size. 
Nevertheless, there could be several copies of the Universe within the horizon. If the 
Milky Way was once a quasar we might detect its 'ghost' images. Using new large 
quasar catalogues we repeat the search by Fagundes & Wichoski for antipodal quasar 
pairs. By applying linear theory to account for the peculiar velocity of the local group, 
we are able to narrow the search radius to 134 arcsec. We find seven candidate an- 
tipodal quasar pairs within this search radius. However, a similar number would be 
expected by chance. We argue that, even with larger quasar catalogues, and more 
accurate values of the cosmological parameters, it is unlikely to be possible to identify 
putative ghost pairs unambiguously, because of the uncertainty of the correction for 
peculiar motion of the Milky Way. 
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quasars: general - cosmology: observations - cosmology: large-scale structure of the 
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1 INTRODUCTION 

The idea that we inhabit a compact multi-connected Uni- 
verse has been presented by ma ny authors (see review by 
iLachieze-Rev fc Lumine3 1^9^). The simplest such Uni- 
verse is one in which space is flat and may be visualised 
as tiled by a 'fundamental parallelepiped' with lengths Ox, 
ay, Oz- If the sides of the fundamental domain are signifi- 
cantly smaller than the horizon length, Kh, observable ef- 
fects would be present. One of the most interesting, albeit 
surreal, effects would be the appearance of 'ghost images'. 
This is where radiation from an object, such as a galaxy, 
has traversed the fundamental domain, resulting in a repeat 
or ghost image of the object. The simplest ghost images to 
search for are ones of our own galaxy, since they appear at 
the nodes of a lattice of which the Milky Way is the origin. 
Therefore ghost images might be most simply identified by 
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searching for antipodal pairs. In this paper we merge three 
new large quasar catalogues to search for ghost images of 
the Milky Way. 

The best limits on the size of the fundamental domain 
come from analyses of the COBE measurement of the cos- 
mic microwave background. The simplest space is the hyper- 
torus. This is a spatially fiat topology where opposite faces 
of the parallelepiped are identified. There are three other 
compact spaces, involving twists of n/2 or tt in pairing op- 
posite faces of a parallelepiped (e.g. lLevinl fl998l. There are 
also two compact models in which the fundamental poly- 
hedron is a hexagonal prism w ith twists of 7r/3 and 2n/3 
iLachieze-Rev fc Luminetll995t) . For equal sides of the par- 
allelepiped~ lLevinri^9^ obtained a minimum size of the 
side of the fundamental dom ain of 0.8Rg, and a similar re- 
sult fo r the hexagonal prism. iBond. Poeosvan fc Souradeenl 
i200(Jl obtained a somewhat larger limit for the hypertorus. 
For unequal sides the constraints are less strong, and there 
remains the possibility that there are many copies of the 
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Universe within the horizon. The existence of antipodal pairs 
is a common feature of most topologies whether the geome- 
try is hyperbolic, flat or spherical. To this extent, therefore, 
although we assume a flat geometry our analysis is also rel- 
evant to nearly-flat geometries. 

On the basis of these COBE limits, it is clear that 
any search for ghosts needs to use objects at large red- 
shifts, and so catalogues of quasars are an obvious source. 
Man y authors have used quasars to probe topological struc- 
ture llFagundes fc Wichoskilll987l iFaeundejllQsllRoukemal 
Il99d^ . but have yet to constrain topology using these meth- 
ods. Following the suggesti on that many large galaxies 
go through a quasar phase, iFagundes fc Wichoskil lll987l) 
searched a quasar catalogue for antipodal quasar pairs in 
the hope of finding ghost images of the Milky Way. In fact 
even with complete catalogues of quasars over the whole 
sky, good fortune would be required for success with such a 
search. This is because quasar lifetimes must be short rela- 
tive to the age of the Universe, so the method relies on the 
quasar being alive at the particular time corresponding to 
the look-back time of any copies. Nevertheless quasar cat- 
alogues have grown enormously in the last few years, due 
in particular to the contributions of the 2-degree-field QSO 
survey (hereafter 2dF) and the Sloan Digital Sky Survey 
(hereafter SDSS). Therefore we considered it timely to re- 
peat the search of Fagundes & Wichoski for antipodal quasar 
pairs. 

The outline of the paper is as follows. In §2 we describe 
the quasar catalogues used. In §3 we explain the search 
method. The search is complicated by the peculiar velocity 
of the Milky Way Vo , which displaces the ghost images from 
opposition by order Vo/c, and so must be accounted for. A 
suitable search radius, then, depends on the uncertainty of 
the historical motion of the Milky Way, which depends, in 
turn, on the uncertainty of the cosmological density param- 
eter. The results of the search are presented in §4, together 
with a brief discussion. We assume flat Universal geometry 
(ilm + f^A = 1) throughout and consider only topologies 
without axial twists. 



2 QUASAR CATALOGUES 

The quasar catalogue used for this investigation was com- 
piled from three sources: 

(i) The complete candidate list for the 2dF QSO survey 
(Smith et al., 2003, in prep.). This catalogue contains 47,768 
objects, from which we extracted the 23,340 quasars with 
spectroscopic redshift s. A subset of 10,00 of the quasars 
has been published in ICroom et al.l ||20^) . 

(ii) The early data rel ease from the Sloan Digital Sky 
Survey Quasar Catalogue dSchneider et alj200ll) which con- 
tains 3,814 quasars. 

(iii) The quasar catalogue of IVeron-Cettv &: VeronI 
which contains 23,760 quasars. 

Duplicat e qua sars, primarily the 10,000 quasars in 
ICroom et all ll200j) which are included in the Veron cat- 
alogue, need to be eliminated. We searched the merged list 
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Figure 1. Sky distribution of the 43,271 quasars in the merged 
catalogue used for the search. 



for pairs of objects within 2 arcsec and Az < 0.02, and in 
each case eliminated one of the pair. This could, in fact, have 
included real pairs, but this would not affect our search for 
antipodal pairs, since it uses a much larger search radius. 
This left a final catalogue of 43,271 objects reaching a maxi- 
mum redshift of 2; = 5.8. The sky distribution of the quasars 
in this merged catalogue is illustrated in Figure 



3 METHOD 

3.1 Peculiar Velocity 

Because of the past motion of the local group, ghost im- 
ages of the Milky Way will appear displaced from the lat- 
tice nodes. This angular displacement is of order Vo/c ^ 7 
arcmin, and must be corrected for. The effect is greatest for 
objects lying perpendicular to the direction of local group 
motion, and falls to zero for objects that lie on the path of 
local group motion. Therefore the idea is to apply an appro- 
priate angular correction for each object in our catalogue to 
bring true antipodal ghost images into opposition. The un- 
certainty of the computed correction sets the search radius, 
since it is much greater than the uncertainty of the quasar 
coordinates, or any other relevant effect. Previous searches 
for antipodal quasar pairs have failed to correct for the mo- 
tion of the local group. 

To find this angular correction we turn to linear theory 
for an expression for peculiar velocity. From this we can 
find the comoving distance moved by the local group as a 
function of redshift, and from this compute the appropriate 
angular correction for each quasar. 

The peculiar velocity, Vpec, is related to the peculiar 
acceleration, g, by 

_ a{z) . _ 2/ 

l|Peeblejll993h where a(z') is the scale factor at epoch z and 
a; is a comoving velocity. Here 

t^ = an(^)"'(l + 2)'. 



available at: http: / /www.obs-hp.fr/www/catalogues/veron2_10/veron2_10.html 



The term / is the dimensionless velocity factor, given by 
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Figure 2. The angular correction to the quasar position as a func- 
tion of redshift, for a quasar located in a direction perpendicular 
to the direction of local group motion (maximum correction) . 
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where 5 is the overdensity or density contrast. For a flat 
universe, a good approximation to / is given by 

jLahav et alJll99ll) . 

In our hnear approximation the pecuhar acceleration is 
given by 

^ _ GMp{Xp - Xlg) 

where a; is a comoving distance and the subscripts refer to a 
single point-attractor, p, and the local group, Ig. This equa- 
tion provides the scaling of g with a, but the mass of the 
attractor, and its distance, are not needed explicitly, since 
we have a measurement of the peculiar velocity today. 

The expression for peculiar velocity is numerically in- 
tegrated to find the comoving distance moved by the local 
group as a function of redshift. This is then converted to an 
angular correction for a ghost image at redshift z. The only 
input variables to this correction algorithm are Jim and Vo- 
The values adopted for this analysis are: 

(i) nm = 0.3±0.1 (e.g.l^rdeeTaOEQoi). 

(ii) Uo = 627 ± 22km s-i i n direction (l b) ^ (276° ± 
3°, 30° ± 3°) from COBE data llKogut et alJIlQOaT) . 

The angular correction for each object depends on red- 
shift and the angle between the quasar, the observer, and 
the direction of local group motion. For an object located 
perpendicular to the direction of local group motion (max- 
imum correction) this correction is approximately 7 arcmin 
at 2 = 0, and increases to approximately 12 arcmin at z = 5. 
A plot of this angular correction is shown in Figure |5| 

3.2 Search Radius 

An appropriate search radius must be related to the un- 
certainties in the input data, and also involves a trade-off 
between contamination and incompleteness. Too large a tol- 
erance and we recover too many false positives, swamping 
any real pairs, while with too small a tolerance we risk losing 
the few putative real pairs in the data. 

We required that the redshifts of the pair agree within 
a tolerance = 0.04. The large majority of the redshifts in 
the catalogue have been measured from intermediate reso- 
lution spectra, many by cross correlation against a template 



spectrum. For these the accuracy of th e redshift is better 
than = 0.01 (e.g. lCroom et al.ll200ll) . Therefore the cho- 
sen tolerance will include nearly all real pairs in the cata- 
logue. 

The position of any quasax in the catalogue, after cor- 
rection for the peculiar velocity of the local group, is uncer- 
tain, because of the uncertainty of the parameters flm and 
Vo- Gravitational lensing could also contribute, since the 
light paths from antipodal quasars traverse different regions 
of the fundamental domain. This is a much smaller source 
of uncertainty, however, and may be neglected. It would be 
very unlikely for a ghost pair to be lensed out of opposition 
by more than a few arcsec. 

In principle we could compute a probability distribu- 
tion function (pdf) for the coordinates of each quasar, and 
select a contour level at, say, the 95 per cent confidence level. 
Then, in searching for an antipodal quasar, the size of the 
contour should be doubled. This is because the error in the 
peculiar velocity correction (due to poor choice of Sim, Uo) 
is identical for each quasar of an antipodal pair. This proce- 
dure would be computationally expensive, and we followed 
a simpler Monte Carlo approach. The algorithm to correct 
the coordinates was run on the catalogue 200 times, each 
time using values of the input parameters fim, Uo, chosen 
from a Gaussian distribution about their mean. We com- 
puted the angular displacement from the mean position for 
each quasar in every simulation, and selected the angle which 
included 95 per cent of all the points. This angle was 67 arc- 
sec. Therefore the procedure for the search was to correct 
the quasar coordinates using the best estimates for 57ni, and 
Vo, and then for each quasar to search for antipodal quasars 
in the catalogue within a search angular radius of 134 arcsec. 
Because the coordinate pdfs for each quasar differ in size, 
shape, and orientation, in some cases this search radius will 
be unnecessarily large, and in some cases too small. But, by 
construction, it will include 95 per cent of any real antipodal 
pairs. 

We found that the largest source of error is from the 
uncertainties in our current knowledge of the speed and di- 
rection of the local group relative to the CMB. Our motion 
within the local group is a source of error unaccounted for in 
our search radius determination. As we are searching for a 
ghost of the Galaxy rather than a ghost of the local group, 
our motion within it should be considered. However, this 
motion is small compared with the large scale streaming 
motion of the local group itself. 



4 RESULTS AND DISCUSSION 

With a search radius of 134 arcsec and Az = 0.04 we re- 
cover seven antipodal quasars pairs. Details of the pairs are 
provided in Table The table includes the coordinates of 
each quasar in the pair both before and after correction for 
peculiar velocity as well as the redshifts, brightnesses, the 
catalogue containing the quasar, the redshift difference, and 
the angular displacement from opposition of each pair. 

The relative positions of the two quasars in each pair 
are plotted in Figure For each pair the plot is centred on 
the antipode of the first quasar in the pair, and the points 
illustrate the offset of the second quasar from perfect anti- 
alignment, for the 200 Monte Carlo realisations of the pa- 



4 S. J. Weatherley, et al. 



.*Q 0007+0236 
200 2QZJ121050.8-G24456 



100 



-200 - 



100 -100 -200 
Aa (orcsecs) 



100 - 



-100 - 



-200 



*QNY1;1 1 




003551.98 + 005726.3 

























200 100 -100 -200 
Aa (orcsecs) 



200 



100 



-200 



*Q 0058-3424 
*KP 1259.4+34.7 


• , • . * • 

. J* • -;. • 







200 100 -100 -200 
Aa (orcsecs) 



2240.0+0019 
2QZJ104229. 5-001616 . 



100 



-200 



200 



100 



-100 



-200 



• ..■■All*- ^ 
•'•",*-tf 



200 



100 -100 -200 
Aa (orcsecs) 



2QZJ1 12525.9-010342 
232519.33 + 01 1 147.8 




t ' 





200 



100 -100 -200 
Aa (orcsecs) 



-FIRST J1057 + 31 19 
200 2QZJ225706. 5-310312 



100 



-100 



-200 - 




200 100 



-100 -200 



Aa (orcsecs) 



100 - 



-100 - 



-200 



2QZJ133904. 3-001945 
013837.29 + 002818.5 









200 100 -100 -200 
Aa (orcsecs) 



Figure 3. Plots of the seven candidate ghost pairs indicating the degree of anti-alignment, and its variation with the input parameters. 
For each pair the cross is centred on the antipode of the first quasar, and the triangle marks the displacement of the second quasar 
for the best values of the parameters Qui, and Vo- In each case the second quasar lies within the search radius of 134 arcsec, which is 
how they were selected. The dots show the computed displacements for the Monte Carlo variations of the parameters Qm , and Vo ■ This 
distribution overlaps the antipode in each case, indicating that for possible values of the parameters Qm , and Do the two quasars could 
be perfectly anti-aligned. 



rameters Qm, and Vo- In each case the distribution of points 
overlaps perfect anti-ahgnment, indicating that, given the 
uncertainties in the parameters, the pairs could be real ghost 
images. As a further investigation, the angles between the 
axes of each quasar pair were computed to search for any 
specific angles that could indicate points on a lattice (such as 
7r/2), but nothing striking was discovered. There are magni- 
tude differences, up to Am = 2, between the two quasars of 
a pair, but these might be due to variability because of the 
epoch differences of the observations, or time delay differ- 
ences because the light paths to us traverse different regions 
of space, or the differences could be because each quasar is 
being viewed from a different direction. For the same rea- 
sons it would not be unreasonable for the spectra to differ 



also. A more important consideration is the number of pairs 
that might be expected by chance. 

We investigated how the expected number of antipodal 
pairs varied with respect to the search tolerance. We varied 
the search angular radius from G = to = 10 arcmin, and 
the redsliift tolerance from Az = to Az = 0.4. Although 
the catalogue is inhomogeneous we found that the number 
of antipodal quasar pairs recovered followed closely the rela- 
tions DC 0^ and oc Az expected for unrelated pairs. Therefore 
we used a fit to the relations at large Q, Az to establish an 
estimate of the number of pairs expected by chance. For our 
search tolerance of = 134 arcsec we would have expected 
12 ± 3.5 quasar pairs by chance. Our result of seven pairs, 
therefore, is consistent with chance. While it remains possi- 
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Table 1. Antipodal pairs of quasars within 134 arcsec and Az = 0.04. Redshifts preceded by J are of lesser accuracy or even wrong (see 
Vc5ron catalogue for further details). 



Name 


Position (J2000) 
R.A. 


Dec 




Corrected Position 
R.A. 


Dec 




Redshift 


Magnitude 


Catalogue 


Deviation from 
opposition 


*Q 0007+0236 
2QZJ121030. 8-024456 


OOli 10m 19.5s 
12h 10m 30.88s 


+02° 53' 
-02°44' 


'38" 
'56.3" 


OOh 10m 23.1s 
12h 10m 27.2s 


+02° 49' 
-02°48' 


57.5" 
'39.7" 


0.588 
0.6257 


V— 18.4 
bj=20.1 


^Aeron 
2dF 


6 — 99" 

Az = 0.0377 


*Q 0058 — 3424 
*KP 1259.4+34.7 


Olh 01m 14.1s 
131i 01m 46.5s 


-34°08' 
+34°26' 


'31" 
'34" 


Ollr 01m 34.7s 
13h 01m 25.8s 


-34° 16' 
+34° 18' 


'58.0" 
06.0" 


^2.08 
t2.08 


V=19.36 
V=18.5 


\^eron 
Vcron 


9 — 129" 
Az = 0.0000 


FIRST J1057+3119 
2QZJ225706. 5-310312 


lOh 57m 5.2s 
221i 57m 6.56s 


+31°19' 
-31°03' 


'07" 
'12.5" 


lOh 57m 06.7s 
22h 57m 05.0s 


+31°ll' 
-31°11' 


01.5" 
'l7.2" 


1.332 
1.3398 


V=18.7 
bj=20.1 


Vcron 
2dF 


B = 27" 

Az = 0.0078 


*QNY1:11 
003551.98+005726.3 


121i 36m 9.9s 
OOh 35m 51.98s 


-00°48' 
+00°57' 


'03" 
'26.3" 


12h 36m 02.8s 
OOh 35m 59.0s 


-00°52' 
+00° 52' 


'56.0" 
33.9" 


1.874 
1.905 


V=21.3 
g/=19.06 


Vcron 
SDSS 


= 60" 

Az = 0.0310 


*2240. 0+0019 
2QZJ104229. 5-001616 


221i 42m 33.3s 
lOli 42m 29.53s 


+00°27'09" 
-00°16'l6.l" 


22h 42m 31.6s 
lOh 42m 31.3s 


+00° 22' 
-00°20' 


30.3" 
'55.7" 


t2.0 
1.9713 


V=20.0 
bj=19.7 


Vcron 
2dF 


e = 95" 

Az = 0.0287 


2QZJ133904. 3-001945 
013837.29+002818.5 


131i 39m 4.35s 
Olh 38m 37.29s 


-00° 19' 
+00°28' 


'45" 
'18.5" 


13h 38m 50.6s 
Olh 3Sm 50.9s 


-00°24'44.3" 
+00°23'20.2" 


1.3474 
1.35 


bj=20.4 
g/=18.93 


2dF 
SDSS 


e = 84" 

Az = 0.0026 


2QZJ112525. 9-010342 
232519.33+011147.8 


llh 25m 25.95s 
231i 25m 19.33s 


-01°03' 
+01°11' 


'42.4" 
'47.8" 


llh 25m 24.7s 
23h 25m 20.6s 


-01°08' 
+01°07' 


'07.6" 
24.4" 


1.743 
1.724 


bj=20.4 
g/=18.45 


2dF 
SDSS 


e = 74" 

Az = 0.0190 



ble that one or more of the pairs could be a real ghost pair, 
this cannot be proven from the existing data. 

4.1 Discussion 

We have presented a search for ghost images of our galaxy 
using the method of antipodal quasar pair detection, correct- 
ing for the angular distortion due to the history of motion 
of the local group. We have found no significant difference 
in the number of pairs found over the number we would 
expect purely by chance using our search tolerance on this 
catalogue. We have examined the angles between the axis of 
each quasar pair but failed to discover any indication of a 
lattice structure. 

The uncertainties both with our present knowledge of 
cosmological parameters and the peculiar velocity of the lo- 
cal group make a search such as this increasingly difficult as 
quasar catalogues increase in size. To prove significance in a 
search such as this the search tolerance needs to be reduced 
to a value where we would be expecting only one or two de- 
tections purely by chance. In the future we can expect im- 
provements in the accuracy of the parameters Sim, Vo, that 
will allow a smaller search radius, that could compensate 
for the much larger future quasar catalogues. However, soon 
the uncertain correction for the peculiar motion of the Milky 
Way within the local group will become a significant factor. 
Given the COBE constraints on the size of the fundamental 
domain (and so the small possible number of ghost images) , 
and the relatively short lifetime of quasars, it would seem 
unlikely that it will be possible to identify putative ghost 
images, unambiguously, by this method. Instead, extension 
of the search to configurations of three or more quasars on 
a lattice could yield a convincing detection. 
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